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We study the ± J three-dimensional Ising model with a spatially uniaxially anisotropic bond ran- 
domness on the simple cubic lattice. The ± J random exchange is applied in the xy planes, whereas 
in the z direction only a ferromagnetic exchange is used. After sketching the phase diagram and 
comparing it with the corresponding isotropic case, the system is studied, at the ferromagnetic- 
paramagnetic transition line, using parallel tempering and a convenient concentration of antifer- 
romagnetic bonds (p z — 0;p xy — 0.176). The numerical data point out clearly to a second-order 
ferromagnetic-paramagnetic phase transition belonging in the same universality class with the 3d 
random Ising model. The smooth finite-size behavior of the effective exponents describing the 
peaks of the logarithmic derivatives of the order parameter provides an accurate estimate of the 
critical exponent 1/v = 1.463(3) and a collapse analysis of magnetization data gives an estimate 
/3/v = 0.516(7). These results, are in agreement with previous studies and in particular with those 
of the isotropic ± J three-dimensional Ising at the ferromagnetic-paramagnetic transition line, indi- 
cating the irrelevance of the introduced anisotropy. 

PACS numbers: 75.10.Nr, 05.50.+q, 64.60.Cn, 75.10.Hk 



I. INTRODUCTION 

Ising spin glass models yield phase diagrams with dis- 
tinctively complex ordered phases, in d = 3. These mod- 
els, although relatively simple in their formulation, have 
been proposed to describe complex systems exhibiting 
frustration, e.g., materials such as Fei_xiMn X jTi03 and 
Eui_ x Ba x Mn0 3 [Ml, neural networks [J] etc. 

The simplest of such models, but most influential over 
the years, is the Edwards Anderson model [5|, Q defined 
by the Hamiltonian 

H = - ^JjjSjSj, (1) 

m 

where the summation is over nearest-neighbors, s = ±1 
and Jij denotes the uncorrelated quenched exchange in- 
teraction. There are two popular quenched random dis- 
order distributions, the Gaussian distribution of random 
bonds with zero mean and unity standard deviation and 
the bimodal distribution of Jij given by 

P(J l3 ) = pd(J t3 + 1) + (1 - p)8(J ij - 1) (2) 

Recently, the spatially uniaxially anisotropic d = 3 
spin glass system has been solved exactly on a hierar- 
chical lattice by Guven et al. Q. Their general study 
revealed a rich phase diagram topology and several new 
interesting features. The Hamiltonian of this anisotropic 
variant, differentiates, along z axis, the probability dis- 
tribution of quenched randomness, but also the strength 
of the exchange interaction and can be written 



Accordingly, the bimodal distribution of J," takes the 
more general form 

P(Jg.) = Pu 5(J^ + J u ) + (1 - p u )6(J% - J%) (4) 

where u denotes the z axis (u = z) or the xy planes (u — 
xy), J," denote the corresponding exchange interactions 
and p u are the probabilities of two neighboring spins (ij) 
having antiferromagnetic interaction along z axis or the 
xy planes. 

The standard isotropic case, defined by Eq. ([1]) and 
Eq. ©, corresponds to Jfj = J* v and p z = p xy . The 
global phase diagram of this isotropic case, is sketched 
in Fig. [T]and, as shown, there exist (in d = 3) three dis- 
tinctive phases, ferromagnetic, paramagnetic and glassy 
phases. All transitions among these phases, are be- 
lieved to be of second order and also to belong to dif- 
ferent universality classes. Several accurate studies have 
been carried out, to determine the critical behaviors 
along these transition lines, for the finite-temperature 
phase transitions [8l4l7|. Most of these studies concern 
the ferromagnetic-paramagnetic (FP) and the glassy- 
paramagnetic (GP) lines. There is also a recent study 
for the ferromagnetic-glassy transition line Q. The 
transition lines meet in a multicritical point [18l423l | lo- 
cated along the Nishimory line [23-|27l w ith coordinates 
T M = 1.6692(3), p M = 0.23180(4) \m$2Aj. The FP tran- 
sition line, starts at the pure Ising model (p — 0), for 
which all critical properties have been extensively stud- 
ied and a recent very accurate estimate of the correla- 
tion length exponent is v — 0.63002(10) [2l j| -[3 1 f with a 
critical temperature T c = 4.5115232(16) (25(31]. As 
shown by Hasenbusch et al. @, the introduction of the 
±J quenched randomness changes the universality of 
the model to that of the random Ising model (RIM) 
or randomly diluted Ising model (RDIs) [j| [l6|, f34j . in 
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which several spin models appear to belong. These, in- 
clude models such as the randomly site and bond di- 
luted, the random bond [35M39( | . the random bond 
Blume-Capel |40| model and, of course, the already 
mentioned isotropic ±J three-dimensional Ising at the 
ferromagnetic- paramagnetic transition line Q . An accu- 
rate estimate of the correlation length exponent for these 
d = 3 FP transitions, characterizing the RIM universality 
class, is that given by [35j v = 0.6837(53). 

For the glassy-paramagnetic (GP) transition line (MB 
line in FigfJL most of the work has been curried out 
at p = | |ll|, MM, EH-ii. However, because of 
severe inherent difficulties, due to both strong frustra- 
tion and disorder effects in this region, there is a large 
spread in the estimates of critical exponents and there 
remain some questions related to universality. This is 
very clearly reflected in Table I of Katzgraber et al. [IIJ 
in which one can observe a very large spread in the cor- 
relation length exponent, and even in the estimates of 
critical temperature. We shall quote from this paper, 
the estimates v — 2.39(5) and T c = 1.120(4), which ap- 
ply to the present isotropic ±J spin-glass model. We 
will also quote the more recent estimates v = 2.45(15), 
T c = 1.109(10) UM. Finally, from the recent work by 
Ceccarelli et al. Q , we know that the correlation length 
exponent for the ferromagnetic-glassy (FG) transitions 
has been estimated to be v — 0.96(2). 

Several other features of the critical behaviors and 
the global phase diagram are known for the isotropic 
case. Ground state calculations by Hartmann [54| (pa = 
0.222(5)), indicated a reentrant FG transition line. This 
was nicely verified by the finite temperature study of 
Ceccarelli et al. since they estimated p(T = 0.5) = 
0.2271(2), predicting the ferromagnetic-glassy transition 
line to be slightly reentrant. These results are in ac- 
cordance with the Nishimori expectations [25J, |26| that, 
this line cannot be forward and are also reflected in 
renormalization-group calculations @, [H, [H| . 

In this paper we will focus on the nature of the FP 
transition for the anisotropic model described by Eq. ([3]) 
and Eq. (TJ| by considering the particular model case p z — 
0',Pxy < 0.5 with Jfj = J^ v . The present study is part of 
a research program to study by Monte Carlo methods the 
general spatially uniaxially anisotropic d = 3 spin glass 
system considered by Guven et al. [7| . The main motiva- 
tion is to provide numerical evidence for the universality 
question, by investigating possible effects caused, by the 
introduced anisotropy, on the critical exponents along 
the three different transition lines. We hope that our 
attempts may also play some role in future studies and 
will provide new information for the critical behavior of 
frustrated systems. 

We start by describing our current perception on the 
global phase diagram of the above anisotropic case. This 
phase diagram is sketched also in Fig.[TJ by the dash lines. 
As shown, there exist again(in d = 3) three distinctive 
phases, ferromagnetic, paramagnetic and glassy phases. 
The points shown on the FP transition line have been 




FIG. 1. (color on line) Phase diagram points (full squares, 
asterisk, and dot) of the isotropic Edwards- Anderson model 
are taken from |j, [Til . ITU , [54| . The phase diagram is drawn 
in solid lines and separates the three phases: ferromagnetic 
(F), spin-glass (G), and paramagnetic (P). These lines meet 
at the multicritical full dot point (M). Dash lines (and open 
symbols) illustrate the phase diagram of the anisotropic model 
p z — 0,p — p X y The open triangles, asterisk and dot points 
are discussed in more detail in the text. In particular, the full 
asterisk (p = 0.117) corresponds to the improved model case 
of [9( and the open asterisk (p = p xy — 0.176) to the improved 
case studied here. 



located from the peaks of the magnetic sample-averaged 
susceptibility for a lattice of size L = 16. The open as- 
terisk (see also the solid drop lines) on this line indicate a 
special point corresponding to p* xy = 0.176. This choice 
was inspired by the improved model case analyzed by 
Hasenbusch et al. Q for the isotropic model (see on the 
isotopic phase diagram the case with p* = 0.117 and the 
corresponding dash drop lines). According, to Hasen- 
busch et al. [9J for the isotropic case at p* — 0.117(3) 
the leading scaling corrections vanish and this case pro- 
vides an improved model for accurate exponent estima- 
tion. Our quest for a convenient and improved model case 
for the anisotropic FP transition, directed as to find the 
phase diagram point corresponding to about the same 
temperature as the isotropic phase diagram point. As 
can be seen, by comparing the drop lines, this would 
suggest the case p* xy = 0.176. Alternatively, let us try to 
satisfy the relation (1 - p* y (T)) = 1.5(1 - p*{T)) - 0.5 
between the two FP transition lines. This prognostic 
method assumes that phase diagram points correspond- 
ing to the same temperature would have (approximately) 
the same ratio of ferromagnetic to antiferromagnetic in- 
teractions. Putting p* — 0.117(3) in this relation, we 
find p* xy = 0.1755(30), which is practically the above 
value for the anisotropic case. As it will be seen in 



3 



the sequel, the proposed case p* xy — 0.176 provides in- 
deed an excellent case and gives, for the anisotropic 
model, a very accurate estimate of the correlation length 
exponent^ = 0.6835(25)). 

The rest of the anisotropic phase diagram (see also the 
caption of Fig. [1}, was mainly constructed on the premise 
that the above relation may hold approximately for the 
lines separating the ferromagnetic phase from the other 
two phases. Thus, applying this prognostic to the point 
p(T = 0) = pa = 0.222(5) of the isotropic case, we find 
that p xy (T = 0) = PA ani = 0.333(5). Quite remarkably 
we found from our ground state calculations a value that 
perfectly verify this prediction. Based on the rather short 
calculation of ground states on lattices sizes L = 6, 8, 10 
we found, by a collapse method, PA ani — 0.332(12) aston- 
ishingly close to the predicted value. After this, the mul- 
ticritical point for the anisotropic case, shown in Fig. [T] 
has been placed by a mere application of the above re- 
lation giving pM ani = 0.34770(4). One should, of course, 
realize that this prognostic is only an approximation, 
which could not be expected to apply in the more gen- 
eral case. It is possible that, it is only successful for 
the present case giving a reasonable approximation for 
the FP and FG transition lines. To complete our per- 
ception of the anisotropic phase diagram, we point out 
that for the case p z = 0;p X y < 0.5 we found from Monte 
Carlo simulations a phase diagram point which appears 
to coincide with the isotopic phase diagram point B. For 
instance, extrapolating the L = 6,8, 10 Binder's fourth 
order cumulant crossings of the overlap order parame- 
ter we found Tg arii = 1.07(4), very close to the cor- 
responding estimate of the isotropic case. This finding 
appears to be a very interesting and physically appeal- 
ing prediction. A proof that Tg jso = Tg ani exactly, 
may also help a better understanding of the general uni- 
versality question. In closing this discussion, we shall 
mention that the introduction of anisotropy changes in 
general the phase diagram and possibly the symmetry 
of multicritical points, as discussed by Guven et al. [7[. 
These multicritical points, such as in the present model 
case (p z = 0;p xy < 0.5), do not necessarily obey Nishi- 
mori conditions and one can not exclude the possibility 
of a ferromagnetic-glassy transition line being forward. 
However, our prognostic method, insist that the present 
anisotropic FG transition line is reentrant. 

The rest of the paper is laid out as follows: In the 
following subsection we give a description of our nu- 
merical approach utilized to derive numerical data for 
large ensembles of realizations of the disorder distri- 
bution and lattices with linear sizes within the range 
L € {8 — 44}. Then in subsection lll Bl the finite-size scal- 
ing (FSS) scheme is described in some detail. In Sec. [Hi] 
we present all our FSS attempts, that give good esti- 
mates of all critical exponents and verify that the present 
anisotropic model belongs to the universality class of 
RIM. Our conclusions are summarized in Sec. IIV1 



II. MONTE CARLO SIMULATIONS AND 
FINITE-SIZE SCALING SCHEME 

A. MONTE CARLO METHOD 

In the present paper, we shall use our recent approach 
to disordered systems [40| , based on a parallel tempering 
(PT) practice. Our PT protocol will use an adequate 
number of Metropolis [57j sweeps of the lattice, as an el- 
ementary Monte Carlo step, so that the correlation times 
of the PT protocol will be, for all temperatures used, very 
close to 1.0. We mention here, that parallel tempering, 
combined with Metropolis algorithm has been used also 
by Ceccarelli et al. [8j in their study of the FG transition. 
Further, as pointed out by Hasenbusch et al. [H , in their 
study of the isotropic case, the Metropolis algorithm may 
be more effective than cluster dynamics, for intermediate 
lattice sizes, as a result of frustration effects present in 
the ± J models. 

Our PT protocol (also our FSS approach) has been de- 
scribed in detail in our recent paper pol [58j j and only a 
brief summary is given here. For the estimation of the 
critical properties we generate MC data that cover several 
finite-size anomalies of the finite systems of linear size L. 
The PT approach is carried over to a certain temperature 
range depending on the lattice size. These temperatures 
are selected in such a way that the exchange rate is 0.5, 
using a practice similar to that suggested in Ref. [59| . 
The appropriate temperature sequences were generated 
via short preliminary runs in which we apply a simple 
histogram method [60h62| | to determine from the energy 
probability density functions the temperatures, satisfying 
the above exchange condition [40l I59 ] . The preliminary 
runs cover several disorder realizations and the average 
over the temperatures sequences provides us with a pro- 
tocol with very small variation in the exchange rate con- 
dition, as one moves from one realization to the other. 

The proposed MC scheme was carefully tested for all 
lattice sizes before its implementation for the generation 
of MC data. The tests included the estimation of the MC 
times necessary for equilibration and thermal averaging 
process applied to a particular disorder realization, and 
involved also the observation of running sample-averages 
of several thermodynamic quantities, such as the mag- 
netic susceptibility, at a temperature close to criticality. 
Provided we use reasonably long MC times for equilibra- 
tion and thermal averaging we find an effective cancella- 
tion of statistical errors from the sample-averaging pro- 
cess, since the sample-to-sample fluctuations are larger 
(by an order of magnitude) than the usual statistical er- 
rors (in the thermal averaging process). 

The present simulation task is, of course, quite consid- 
erable, since we have to satisfy both a good equilibration 
of the system and also sum over a very large number of 
disorder realizations. The number of realizations influ- 
ence the accuracy of our data, and the proper selection of 
temperatures influence their suitability for the locations 
of the finite-size anomalies. Since we are implementing 
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a PT approach, based on temperatures corresponding to 
an exchange rate 0.5, we are selecting proper temperature 
sequences consisting of a number of (say 3 or 5) different 
temperatures and averaging for this set of temperatures 
over a relatively large number (~ 1000) of disorder real- 
izations. This yields a number of (3 or 5) points of the av- 
eraged curves [Z] av , where Z denote the thermal average 
of some thermodynamic quantity, such as the magnetic 
susceptibility. The PT protocol is repeated several times 
(depending on the linear size L) by using new sets of tem- 
peratures (translated with respect to the previous sets) 
and a final dense set of points is obtained. The above 
is an efficient and most importantly quite accurate prac- 
tice, since we finally obtain averaged curves correspond- 
ing to a very large number of realizations, describing with 
high accuracy all the averaged finite-size anomalies of the 
system. Monte carlo data were collected for systems of 
linear sizes L = 8, 12, 16, 20, 24, 28, 32, 36, 40, 44 and the 
FSS analysis was performed to the averaged quantities 
obtained from these MC data. 



B. FINITE-SIZE SCALING SCHEME 

In the standard approach of FSS for a random system, 
a large number of disorder realizations has to be used in 
the summations in order to obtain good sample-averages 
of the basic thermodynamic quantities Z, which are the 
usual thermal averages of a single disorder realization. 
From the disorder averages [Z] av we obtain their finite- 
size anomalies, denoted here as [Z)l_ v . These finite-size 
anomalies will be used in our FSS attempts, following a 
quite common practice [63[. Their temperature locations, 
denoted by T^ z \* will be used in the sequel in our FSS at- 
tempts. Thus, our study concerns the critical exponents 
describing the disorder-averaged behavior and we do not 
attempt a FSS analysis based on sample dependent pseu- 
docritical temperatures. The later is, a more demanding 
alternative approach [64]-[66j], which considers the indi- 
vidual sample dependent maxima (anomalies) and the 
corresponding sample dependent pscudocritical temper- 
atures. Note that, for disordered systems one could make 
in principle a clear distinction between typical and aver- 
aged exponents [H, [66| . 

From the MC data, several pseudocritical temper- 
atures may be estimated, corresponding to finite-size 
anomalies and these are expected to follow a power-law 
shift behavior T [z] , v =T c + b z - The traditionally 

used specific heat and magnetic susceptibility peaks, as 
well as the peaks corresponding to the following logarith- 
mic derivatives of the powers n — 1, 2, and n = 4 of the 
order parameter with respect to the inverse temperature 
K = X/T 



parameter derivative 
d(\M\) 



dK 



(\M\H)-(\M\)(H), 



(6) 



dln(M n ) (M n H) 



(5) 



dK (M n ) 
and the peak corresponding to the absolute order- 



will be located and used in our fitting attempts. 

The behavior of the maxima of the logarithmic deriva- 
tives of the powers n = 1, 2, and n — 4 of the order pa- 
rameter with respect to the inverse temperature, which 
as is well known scale as ~ L 1 ^ with the system size [67j . 
will be seen to provide a smooth root for the estimation 
of the correlation length exponent v. Once the expo- 
nent v is well estimated, the behavior of the values of 
the peaks corresponding to the absolute order-parameter 
derivative, which scale as ~ J_y—P)l v with the system 
size (6?| , gives one route for the estimation of the mag- 
netic exponent ratio /3/v. 

For the estimation of the critical temperature, we shall 
use mainly a simultaneous fitting approach of the several 
pseudocritical temperatures mentioned above. However, 
from the MC data for the disorder averaged magnetiza- 
tion, it is possible to follow an optimum collapse method 
which will provide simultaneously estimates for the criti- 
cal exponents /3/v, X/v as well as the critical temperature 
T c . From the scaling hypothesis 

[<|M|>] ou = M(T,L) = L-^f[(T-T c )L^ v ], (7) 

and the disorder averaged magnetization data, we at- 
tempt in the next section the estimation of critical be- 
havior by using a recently published collapse method that 
makes use of the downhill simplex algorithm 68] . 



III. ESTIMATION OF CRITICAL BEHAVIOR 

As pointed out earlier, the present study uses the finite- 
size anomalies of the sample average of the logarithmic 
derivatives of the powers (n = 1,2,4) of the order pa- 
rameter with respect to the inverse temperature, for the 
estimation of the correlation length exponent v. It will 
be seen that, the finite size scaling behavior of the corre- 
sponding peaks provides for the present model an attrac- 
tive and smooth approach to this estimation. Assuming 
that these finite-size anomalies ([Z]*„) of the disorder 
averages [Z] av , where Z is the thermal average given by 
Eq. ((SJ), scale as ~ L x l v with the system size 67 1 we es- 
timate in figures [2] and [3] this exponent. A very good 
scaling behavior is observed already by using the whole 
size range of our Monte carlo data L — [8 — 44]. The 
simultaneous fitting attempt to the expected power-law 
behavior, illustrated in Fig. [3J gives the estimate shown 
in the panel: X/v = 1.468(6). 

By varying the L m i n , in these simultaneous fitting at- 
tempts, we obtain a sequence of effective exponents de- 
pending on the minimum size used. The behavior of these 
effective exponents is shown in Fig. [3J The smooth, al- 
most perfectly linear, behavior of these effective expo- 
nents enables us to estimate confinedly X/v = 1.463(3). 
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FIG. 2. (color online) FSS behavior of the peaks of the log- 
arithmic derivatives of the powers n = 1, 2, and n = 4 of 
the order parameter with respect to the inverse temperature. 
The estimate for the exponent 1/u = 1.468(6) (shown in the 
panel) is obtained by applying a simultaneous fitting attempt 
to a simple power law in the whole size range L = 8 — 44. 



FIG. 3. (color online) Illustration of the behavior of the effec- 
tive exponents (l/v) e ff. The solid line drawn in the panel, to- 
gether with the dot lines, indicate the critical exponent range 
for the present model. The linear behavior illustrated give 
the accurate estimation 1/u = 1.463(3). For comparison the 
analogous narrow range 1/u — 1.587(1) 28], for the pure 3d 
Ising model is shown in the upper part of the panel. 



The error range of this estimation is indicated in Fig. [3] 
by the dot lines and is compared with the correspond- 
ing 3d pure Ising model for which an extremely accurate 
estimation is available [28} . The present estimate for the 
correlation length exponent v = 0.6835(25), compares 
well with the estimate v = 0.683(3) of Hasenbusch et 
al. [9], for the corresponding isotropic ±J Ising model at 
the ferromagnetic-paramagnetic transition line and is in 
excellent agreement with the estimate v = 0.6837(53), of 
the extensive numerical investigations of Ballesteros et 
al. [35[ for the site-diluted Ising model. The good behav- 
ior of the effective exponents in Fig. [31 giving improved 
error bounds, suggest the present model as a suitable 
candidate for further attempts and refinements of this 
exponent, possibly followingthe more sophisticated ap- 
proach of Hasenbusch et al. Q to extend the Monte Carlo 
data to larger lattice sizes. 

We proceed to calculate the critical exponent ratio 
j/u from the peaks of the sample-averaged susceptibility 
(Mau)- We assume that these finite-size anomalies obeys 
a simple power law: [x\* av — b-L 1 '", and follow again the 
practice of observing the behavior of effective exponents 
by varying the L m i n of the fitting range. The resulting 
sequence of effective exponents (L m in = [8 — 24] ) is illus- 
trated in Fig. SJ As seen from this figure, the behavior of 
these estimates is not linear. Thus, the attempted linear 
fit, appears as an overestimation, but gives clearly the 
underestimated estimate 7/V = 1.9614(28), appreciably 
smaller than that of Ballesteros et al. [35] j/u — 1.963(5) 
and even more than those of Hasenbusch et al. 1.964(1) 
[Ifl for the site-diluted 3d Ising model and 1.964(2) Q for 
the corresponding to the present model isotropic ± J Ising 
model at the ferromagnetic-paramagnetic transition line. 



Therefore, the estimation via the susceptibility peaks is 
here unsatisfactory. 
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FIG. 4. (color online) Illustration of the behavior of the ef- 
fective exponents y/u. The solid line drawn in the panel, 
together with the dot lines, indicate the unsatisfactorily esti- 
mated critical exponent range f/u = 1.9614(28). 

We may now attempt the estimation of the expo- 
nent ratio /3/v, via the scaling behavior of the peaks 
corresponding to the absolute order-parameter deriva- 
tive which is expected to scale as \d(\M\) / dK]* av — 
b ■ Ly i ~^l v . The corresponding effective exponent esti- 
mates, illustrated in Fig. [51 show again a behavior which, 
as seen, is not linear. The illustrated in the figure lin- 
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ear fit gives (1 — f$)/v = 0.9572(54), which by using our 
estimate 1/v = 1.463(3), produces for (3/v the range 
0.5058(84). However, this value is noticeably smaller 
than the value (3/v = 0.518, expected from hyperscal- 
ing and accepting an estimate for 7/1/ of the order of the 
above mentioned literature estimates (say for instance: 
7/1/ = 1.964). The situation in this case can be im- 
proved by attempting the linear fit in Fig. [5] only at the 
last three points. The linear fit in these three points 
(L mm = 16,20,24), gives (l-0)/v = 0.948(3) produc- 
ing, as above, /3/v = 0.515(6). Note also, that the second 
order polynomial fit, shown also in Fig. [5] and applied to 
all five points, gives an estimate (1 — [3)/u = 0.943(6) 
producing now fijv = 0.520(9). These values are in good 
accordance with hyperscaling and the literature estimate 

of j/v. 
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FIG. 5. (color online) Illustration of the behavior of the effec- 
tive exponents (1 — /3)/v. The solid and the dash lines drawn 
through the points indicate linear and second-order polyno- 
mial fits, discussed in the text. The critical exponent range 
for the later fit is indicated by the heavy solid and dot lines 
at (l-/3)/u = 0.943(6). 

The critical temperature will be now estimated by a 
simultaneous fitting approach, using several pseudocriti- 
cal temperatures of the sample average of the quantities 
measured (58l. Ir39j. as outlined earlier. The simultaneous 
fitting is attempted to the expected power-law shift be- 
havior T [z] , v =T c + b z - L- l l v for the six pseudocritical 
temperatures mentioned in the previous section. 

We approach this estimation by simultaneous fittings 
in which we are fixing the exponent 1/v to the appar- 
ently accurate estimate 1/v = 1.463. Following our 
earlier practice of using different fitting ranges by vary- 
ing the L m i n of the fitting range, we obtain a sequence 
of estimates illustrated in Fig|6] The linear fit shown 
in the panel gives an estimate T c = 3.2931(12) illus- 
trated with solid and dot lines in this figure, whereas 
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FIG. 6. (color online) Illustration of the behavior of the ef- 
fective pseudocritical temperatures obtained by applying a si- 
multaneous fit on the shift behavior Ti Z ]* = T c + bzL~ and 
fixing the exponents to the accurate estimate 1/v = 1.463. 
The solid dot lines indicate the critical temperature range 
T c = 3.2931(12) obtained by the illustrated linear fit. 



restricting the fit only to the last three points, corre- 
sponding to L m in — 16, 20, 24, gives a higher estimate 
T c = 3.2945(18). We note here that a completely free fit, 
without fixing any parameter, and following the above 
practice gives T c = 3.2934(8) from the linear fit to the 
five points L mm = 8, 12, 16, 20, 24 and T c = 3.2940(16) 
from the linear fit to the three points L m i n = 16, 20, 24. 
Thus, we will suggest that T c = 3.2938(9) 
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FIG. 7. (color online) Illustration of magnetization data col- 
lapse for lattice sizes L — 12 — 44. The method used is 
discussed in detail in the text and gives as shown in the 
panel an optimum collapse at the values T c = 3.2928(7), 
1/v = 1.466(12), and f3/v = 0.516(7). 

Finally, we present the alternative estimation of criti- 
cal behavior by studying the expected scaling law Eq. ([7]) 
for the order parameter data, here the disorder averaged 
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magnetization data. Using the earlier mentioned collapse 
method, we rescale the y axis to ML?'" and the x to 
(T — T^L 1 /" , and we attempt to observe the optimum 
collapse of the magnetization curves taken from differ- 
ent lattice sizes (L = [12 — 44]). We apply the downhill 
simplex algorithm as developed and implemented in ref- 
erence [68J for the estimation of critical properties and 
their error bounds. As shown in the panel of FigjTl the 
optimum collapse gives f3/v = 0.516(7). This value is 
in very good agreement with the expected value as men- 
tioned earlier. In the panel we also give, the resulting 
estimates for the exponent 1/v = 1.466(12) and the esti- 
mate for the critical temperature T c — 3.2928(7). These 
are in fair agreement with our previous findings. How- 
ever, the estimate f3/v = 0.516(7), appears to be very 
satisfactory. 



the reliable estimates (3/v = 0.516(7), by using the col- 
lapse method of reference [68[ and v — 0.6835(25), from 
the smooth behavior of the logarithmic derivatives of the 
order parameter. We have presented also a conjectured 
global phase diagram, providing interesting predictions. 

Currently, we are carrying out further numerical sim- 
ulations. From these, it seems that, the implemented 
here anisotropy (p z = 0; J?- = Jfj*) is also an irrelevant 
parameter for the other transition lines of the phase di- 
agram (the ferromagnetic-spin glass and the spin glass- 
paramagnetic lines). Finally, we are considering the more 
general case ( J?- ^ J*/)- We hope that, we will soon pro- 
vide further confirmation of the discussed in this paper 
predictions and observe and verify the interesting fea- 
tures, of the global phase diagrams, brought out by the 
study of Guven et al. 



IV. CONCLUSIONS 

The present paper pointed out clearly that, the ±J 
three-dimensional Ising model, with spatially uniaxially 
anisotropic bond randomness, give rise to a second-order 
phase transition belonging in the same universality class 
with the 3d random Ising model. The implemented 
anisotropy appears as an irrelevant parameter for the 
ferromagnetic-paramagnetic transition line. We found 
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